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Abstract 

We generalize the approach of Atiyah and Manton for generating Skyrmion 
configurations from instantons. The result is a series whose parameters are 
found directly from the chiral angle equation. The series converge rapidly 
to the exact solution for a class of the Skyrme-like models (including the 
Skyrme model itself) but describe with less accuracy other types of models. 
We describe the procedure and discuss its advantages and limitations. 

PACS number(s): 11.30.Na, 11.30.Rd, 12.15.Cc. 
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I. INTRODUCTION 



As a non-linear theory of pions, the Skyrme model [|TJ provides an approximate descrip- 
tion of hadronic physics in the low-energy limit. In this theory, the nucleon emerges as a 
bound state of the field, or more precisely as a soliton. The original Skyrme Lagrangian was 
of fourth order in field derivatives and it reached a 30% accuracy with respect to physical 
observables. Since then, we proposed along with other authors @-|5| a number of alternative 
models which preserved the form of the original Lagrangian while extending it to higher 
orders. The result was to incorporate effects due to higher-spin mesons and thus to improve 
the fit of most observables. 

All those models have one handicap in common however, they do not admit exact analytic 
solutions. Very few soliton solutions are known as a matter of fact, with the exclusion of the 
one-dimensional sine-Gordon equation , the KdV equation and some other special cases. This 
is a natural consequence of non-linearity, a necessary constraint for a soliton-like solution to 
exist. 

In the absence of analytical solutions, the only alternative to numerical treatment is the 
use of aptly-chosen analytical forms which provide a sometimes sufficient approximation. 
Apart from greatly simplifying calculation of physical quantities, a great deal of information 
can be extracted from such an analytic form. For instance, symmetries and general behavior 
of the solution are much easier to analysis, and the individual features proper to each Skyrme- 
like model also become more explicit. Then, an analytic form proves useful in the stability 
analysis of the soliton, both classical and quantum, and in the calculation of multi-Skyrmion 
interactions. In particular, Igarashi et al. |§ have analyzed the quantum behavior of the 
Skyrme model soliton on the basis of a family of trial functions, taking account of breathing 
motions and spin-isospin rotations whereas Hosaka et al. have examined to the two 
Skyrmion interactions. 

The exact form of most trial functions is largely arbitrary however, with the sole 
constraint that they have the same asymptotic behavior as the numerical solution (i.e., 
F(r) = Ntt — ar + 0(r 3 ) when r — > and F(r) = (3r~ 2 + 0(r~ 4 ) when r — > oo) and that 
they reproduce the numerical results with a degree of accuracy. A few years ago on the 
other hand, Atiyah and Manton brought forward the idea of generating Skyrmion solutions 
from Yang- Mills instantons ||. Using this approach as a starting point, we have developed 
a systematic procedure in order to parameterize any Skyrmion solution. 

We first present the Atiyah-Manton idea and discuss its successes and drawbacks. We 
then introduce our own systematic procedure for parameterizing the solution and apply it 
directly to a few Skyrme-like models in the N — 1 case of a single Skyrmion configuration 
(Section 3), and in the more general case of N > 1 (Section 4). We finally offer some 
comments on the advantages and limitations of our method and on various ways of improving 
it (Section 5) for models where the chiral angle shows non-smooth behavior due to its 
nonlinear nature, . 

II. THE ATIYAH-MANTON IDEA 

The Atiyah-Manton idea is based on topological similarities between Skyrmions and 
instantons. Before introducing it, let us first review the Skyrme model itself. The pion fields 
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are represented by a unitary SU(2) matrix denoted U(x). In its familiar hedgehog form, it 
is expressed as follows: 



U(x) = exp [ir ■ xF(r)] (1) 

where F(r) is called the chiral angle or profile function of the solution. This field configuration 
constitutes a map from physical space R 3 into the group manifold SU(2) and is assumed 
to go to the trivial vacuum for asymptotically large distances. We therefore impose U(r — > 
oo) — >1. From this last condition, one may derive the existence of a topological invariant 
associated with the mapping. The originality of Skyrme's idea was to identify this invariant, 
i.e. the winding number, with the baryon number. 

Introducing the notation L M = Wd^U, the Skyrme Lagrangian for zero pion mass takes 
the form: 

C = -^fTrL^ + ^Tr[L„ L u f (2) 

where the first term coincides with the non-linear sigma model and the second one acts as 
a stabilizer. 

The Yang-Mills instanton follows from a different path. We let A^(x) be a gauge field 
and x = (x, t) be a vector in R 4 . We may then introduce the unitary SU(2) field U(x) in 
terms of the temporal component of A^(x) by the following equation: 

U-\x)d t U{x) = -A t (x) 

This equation is solved through the method of time-ordered products by: 



U (x) = ±T exp 



A t (x)dt 



(3) 



with U (x, — oo) = —1 and U (x, +oo) = U(x). 

The A t component of the Yang-Mills field may furthermore be expressed in terms of the 
scalar superpotential p: 

* = r-v (4) 

2 p 

The so-called 't Hooft instantons of topological number k are then obtained by setting: 

k \2 

where the Xj are poles in R 4 and the Aj are positive constants acting as scales. 

The Atiyah-Manton idea then follows from identifying the field U(x) — U (x, +oo) with 
the Skyrme field, with B=k. In the B=l case of a single soliton centered at the origin 
(Xi = 0), we set: 
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and find by © and (|): 



F(r) 



71 1 



v Vr 2 + X 2 J 

This function has the same asymptotic behavior as the numerical Skyrmion solution: 

7T + ar + 0(r 3 ) , when r — > 

/3r~ 2 + 0(r~ 4 ) , when r — > oo 



(7) 



with a = — 7r/A and /3 = ir\ 2 /2. The A parameter, which gives a measure of the size of 
the soliton, is determined by minimizing the static mass Ms of the soliton. Applying this 
function to the Skyrme model, we find after minimization M$= 104.1 (F n /e), only 0.9% 
higher than the exact numerical result. 

Despite this relatively impressive success however, this approach leaves something to be 
desired. For one thing, the slopes at r — > and r — ► oo both depend on the same parameter 
A, which means that their values cannot be adjusted independently. More precisely, we have: 

3 

@ = 2~- Secondly, we find that the predictive power of this simple function rapidly decreases 
when we consider more elaborate Skyrme-like models such as those we present below. 

Aiming to improve the fit of the exact solution while preserving the general form of the 
Atiyah-Manton ansatz, we consider an extension of the series. The following form turns out 
to simplify the calculations: 



F(r,X) 



71 



1 



oo 
\ " 



A 2 I 



n+l i 
2 J 7"^L 2 J 



„2[I 



Vr 2 + A 2 ^ 



C2n+1" 



n=0 



r 2 + A 



2\n 



(8) 



where [z] is the integer part of z and c = 1. This solution correspond to a general potential 
of the form: 



p = exp 



-hlQ+ J2 d mQ T 



m=l 



where, Q 



r 2 + t 2 
r 2 + t 2 + A 2 



where for d m — 0, one recovers the expression in (^). 

This new function is still characterized by a scale A, as well as by the coefficients c n , 
acting as weights with respect to the corrective terms. As we will show those coefficients are 
not completely free; their value is determined once a choice has been made for the slopes a 
and j3. The scale A is once again determined by minimization of the mass of the soliton. 



III. THE N=l SOLITONS 

It proves useful to express the Lagrangian of all Skyrme-like models in the form (for the 
hedgehog solution): 

oo 

C = J2 h m a m - 1 [3a + m(b - a)] (9) 

771=1 

where a = r~ 2 sin 2 F,b — F' 2 and where we have made the scale change: r — ► s jfr. The 
static energy is then written as: 
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Mi 



or as: 



s = 4vr(— ) / r 2 dr V h m a m - l [3a + m(b - a)} 

6 J m=l 



M s = 4vr(^)/ r 2 dr[3 X (a) + (b-a)x'(a)\ (10) 
e Jo 



where x( x ) = Em=i h m x m and x'( x ) — • Using the same notation, the chiral equation 
becomes: 

// Nrn// n F' -sinFcosF. ... F' rt2 cosF sinFcosF n 
= x '(a)[F" + 2 2 ] + ax"(a)[-2— + F' 2 — — + ]. (11) 



The Skyrme Lagrangian corresponds to the case x{ a ) = Xs( a ) — a + |° 2 - 

More general models should include higher order terms P| . One could consider the model 
for exponential or truncated geometric series: 

Xi(a) = e a -l =a+ia 2 + ^a 3 + ... 

X2M°) = a -\ a2 + l a3 ---jj aM ( 12 ) 

which corresponds to the choice h m <M = — • 

Another interesting model has been proposed by Gustaffson and Riska [|[ which replaces 
the fourth order Skyrme term by a higher order term of order m: 

X3,m(«0 = a + ^a M 
Some other alternatives are due to Jackson et al M. They are: 

Xa{o) = ln(l + a) + l-a 2 

Xs{a) = \\l-e- 2 «]+ l -a + \a 2 



x J a ) =a + 



« 3 



3 + 2a 



a 3 a 4 



X 7 (a) = a + — h 



3 + 4a 1 + 4a 2 

In this work, we shall limit ourselves to these seven models, taking M =20 for the X2,m 
model and M =20 for the X3,m model. 

In order to fix to the coefficients in the parameterization, we use the chiral angle equation 
in the r — > and r — > oo limits. The first step of our procedure is to express the chiral angle 
F(r) as a Taylor expansion, valid for small r: 

F(r) = Nit + a x r + a 2 r 2 + a 3 r 3 + ... (13) 
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and to insert this expression in the chiral equation for each model. Since the coefficients of 
each of the powers of r on the right-hand side of the equation must vanish individually, we 
find in the case of the Skyrme model (i.e. model Xs{ a )) that a n = for n =even and the 
following relations for the lowest n: 



a 3 = 

30 



4 + a? 



a', 



«5 



1400 



1 + a? 
40 + 32al + 11a? + 7a?' 



1 + a? 



2^3 



Thus all «j coefficients are determined once ot\ is fixed. We repeat this operation for r 
using: 



and find f3 n 



F(r) -- 

for n =odd and: 



qf*2i y«3 /^<4 



+ 



10 
.10 



(14) 



oo, 



(15) 



The second step of the procedure is to adjust the parameters of q to these power series 
expansion using both limits: 



F(r -> 0) 



7T 



1 + C 3 ) , (§ + §C 3 - C 5 - C7) ;{ 



-r + 



A 3 



r — ... 



00 



7T 



>2 — C 3 — C 5)A 2 ( — § + §C 3 + |c 5 — C 7 — C 9 )A 4 



Equating the various powers of this expression with the a, and $ coefficients of equations 
(0) and (|IJD, we find: 



a x A /3 2 1 
c 3 = 1, c 5 = — + - - c 3 , c 7 

7T 7tA^ 2 



a 3 A 3 1 3c 3 



Thus all the free parameters of our function are fixed once we have made a choice for 
the scale A and for the slopes ai and As we have pointed out, the value of those three 
parameters is determined by minimization of the soliton mass M$ (eq.([lO|)). 

The calculations are performed up to order n = 10 in (|]). The results we have reached are 
listed in Table |, along with the margin of error with respect to numerical results whenever 
data was available. We notice that the agreement with numerical results varies within a wide 
range, depending on the model considered. The agreement is best when the order of the 
Lagrangian is low (namely for the Skyrme model and the first two Jackson et al. models), 
but for the other models the series does not seem to converge very fast. This observation 
is related to the appearance of a new feature in the solution as m increases. Numerical as 
well as analytical analysis show that the space derivative of the chiral angle F(r) begins to 
develop a cusp near r=l, exhibiting phase transition. 
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The fact that this feature appears in the r ~ A region signals a limitation in our proce- 
dure. The overall shape of our analytic profile function is largely determined by its asymp- 
totic behavior at r — > and r^ oo, and we have no way of carrying a direct fit of the 
solution in-between. When the solution is smooth (i.e., when the order of the Lagrangian 
is low), this causes no difficulty; but when phase transition occurs we cannot expect strict 
agreement between numerical and analytical solutions to hold for all r. 

Next, we consider the Atiyah-Manton approach applied to the N > 1 spherically sym- 
metric solution. Although they are not the lowest energy solutions (stable solitons), the 
approach leads to interesting results since time-ordering in eq. (|3|) is straighforward in these 
cases. 



IV. THE N=2,3,... HEDGEHOG SOLUTIONS 

Proceeding on to the B=2 solution of the Skyrme field equation, we consider the spher- 
ically symmetric case: 

X\ A 2 



(t-Xi) 2 + r 2 (t-T 2 ) 2 + r 2 

The minimal energy is obtained for Ai = A 2 = A and (Tj — T 2 ) large but finite ||. When 
(Ti — T 2 ) — > oo, the solution has a simple form: 

JWr) = 2tt (l + ^=^j = 2*W(r) (16) 

and its energy turns to be larger only by a minute amount (~ 0.5 ppm). Since the introduc- 
tion a parameter for a finite {T\ — T 2 ) seems to bring very little improvement to the solution, 
we shall build our N > 1 solutions on the solution with (T\ — T 2 ) — > oo. 

Applying this function to the Skyrme model and minimizing the mass with respect to A, 
we find M s =310.8 (i^/e), 1.2% higher than the exact result M s =307.1 (F w /e) . Similarly, 
if we set F/v =3 (r) = 3F/v=i(r),we find the minimal mass to be Ms =633.1 (F n /e), 3.0% 
higher than the exact value. We see that the results seem to decrease in accuracy as iV 
increases. 

We generalize our own approach to the N > 1 problem by setting: 

N 

F iV (r,A) = ^F 1 (r,A n ) 

n=l 



We have considered various alternatives for A n . For instance: 

N 

F N (r,\) = Y,Fi(r,\) = NF 1 (r,\) 



71=1 

N 

F N (r, A) = £ F 1 (r, n\) = F 1 (r, A) + ... + F 1 (r, NX) 

n=l 

In the N = 2 and N = 3 cases, the lowest mass Ms is reached as X n = A and A n = nX 
respectively. We show the results in Table || for the N = 2 case. Comparing the numerical 
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and analytical profiles we have obtained for F(r), we have noticed a slight difference. While 
the latter remains smooth for all r, there is a sharper drop-off of the former near F(r) = 7r. 
This disagreement surely accounts for the poorer results we obtain (see Table [TI|). 

V. ADVANTAGES AND LIMITATIONS OF THE APPROACH 

The reliability of our approach seems to depend largely on the model considered. The 
conditions imposed on the c, coefficients insure a good agreement with numerical results 
when r — > and r — > oo. This agreement is conserved when the solution remains smooth for 
all values of r, as in the case of the N — 1 soliton in the Skyrme and Jackson models (i.e., 
models in which the order of the Lagrangian is small). But discrepancies between analytical 
and numerical results begin to appear when we consider higher-order models and when the 
Skyrmion number N exceeds 1. Smaller scale structure then appear in the intermediate 
region (r — > A), and the Atiyah-Manton scheme does not have any special treatment of 
those cases. 

In the multi-Skyrmion problem, the sharper drop-off of the profile function near F = Ntt 
is smoothed out. This may be caused by the fact that our function includes a single scale, 
A. We have tried the alternative of using two different scales by setting F N=2 (r, A) = 
F N=1 (r, A) + Fjv=i(r, 2A), but with little or no improvement in that respect. 

The same type of arguments go for two-phase models. Although the phase transition of 
the solution is reproduced by our analytical function, this feature is not accurately fitted. 
This type of model actually seems to show chaotic behavior near r ~ A. Numerical as well as 
analytical treatment show that when we solve the differential equation for the chiral angle, 
we find we must finely tune the boundary conditions for r — > for a solution to be reached. 
The reason is that any non-trivial solution must obey (r~ 2 sin 2 F)=l to a very good accuracy 
for the range of values < r < tq = 1 almost independently from the rest of the differential 
equation, otherwise terms of the form (r -2 sin 2 F) m go to or oo for m large. Any account 
to parameterize the phase transition must take that into account, which is not the case in 
our procedure. We find that the cm, /3 2 and A parameters must also be finely tuned, and 
that the slopes at and oo do not always correspond to their numerical values. 

Various ways of improving the fit of the solution may be taken into consideration. For 
instance, we have tried using other trial functions, not based on the Atiyah-Manton ansatz. 
The form of those is rather arbitrary however and, despite somewhat improved results, we 
have still met the same basic limitations. Another alternative we considered was to add a few 
more terms to our series in order to account for the coefficients an, (3x2, etc- Unfortunately, 
besides the rapidly increasing complexity of our calculations, we have found that those terms 
contribute little to the r ~ A region. 

Since it seems impossible to carry out a direct fit of the solution in this intermediate 
region, another approach we are considering is to modify the corrective terms of the series, 
attempting to find a pattern in the values taken by the Cj coefficients that would make the 
series converge to a closed form. 

Even though imperfect, we stress again that this approach is still probably the best 
alternative to complete numerical treatment, and that it can prove very useful whenever 
an analytical form of the solution is required, as we have already discussed with regard to 
stability analysis. The approach Igarashi et al. followed with the Skyrme model can easily 
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be generalized to other cases (two-phase structure models, multi-Skyrmion configurations, 
...) by using a new family of trial functions based on our prototype. 
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TABLES 



TABLE I. The values of the a\ and coefficients (see equations 13 and 15) and of some 
physical observables we have evaluated with our analytical ansatz in the N = 1 case (single 
Skyrmion). The Skyrmion mass Ms (see equation 10) and the electrical radius < r 2 >V 2 are 
expressed in units of ^f- and -gp- respectively (with e and F n defined in (2)), the pion decay 
constant in MeV and the axial coupling constant qa is adimensional. We also give the margin of 
error with respect to exact numerical results as found in references [2-5]. 



1 \ 


ai 


O 

P2 


71 AT 

M s 


771 


9A 


, 2 ^ 1/2 


Xs(a) 


2.837 


1.075 


103.129 


64.59 


-0.612 


0.7479 




(0.10%) 


(0.56%) 


/ r\ r\r\r\ A (\~f \ 

(0.0004%) 


(0.03%) 


(0.10%) 


(0.04%) 


Xi (a) 


1.333 


2.8062 


136.505 


72.984 


-0.713243 


1.28496 








^u. 1U /o ) 


^u.u.^ /o ) 


\U.<±0 /o ) 


^U.dU /0 ) 


X2,2o(a) 


0.981 


2.050 


73.5978 


99.29 


-0.743 


1.2090 








(3.59%) 


(3.86%) 


(20.9%) 


(4.27%) 


X3,2o(a) 


0.970 


2.440 


98.3687 


85.42 


-0.766 


1.2778 




(1.25%) 




(0.89%) 


(0.75%) 


(8.27%) 


(1.47%) 


X4(o) 


2.634 


0.982 


84.6796 


71.02 


-0.697 


0.7542 








(0.01%) 






(0.03%) 


X5(a) 


2.587 


1.034 


88.2485 


70.54 


-0.703 


0.7784 








(0.02%) 






(0.24%) 


Xe(a) 


2.933 


0.985 


86.1798 


70.59 


-0.697 


0.7634 








(0.40%) 






(0.07%) 


X7(a) 


3.023 


1.062 


96.6275 


67.55 


-0.663 


0.7824 








(0.65%) 






(0.18%) 
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TABLE II. The values of the ol\ and coefficients (equation 13 and 15), of the Skyrmion mass 
Ms (equation 10) and of the observables F n , g& and < r 2 > 1 / 2 as obtained with our analytical 
ansatz in the N = 2 case. The units are the same as in Table I . The margin of error is given with 
respect to the numerical results we have obtained. 

_ _ _ < r 2 >l/2 

1.472 
2.993 
2.5240 
2.8388 
1.4205 
1.4515 
1.5209 
1.5346 



X{a) 


Oil 


ft 


M S 




9A 


Xs(a) 


5.800 


3.245 


308.569 


40.94 


-0.525 




(4.30%) 


(0.37%) 


(0.48%) 






Xi (a) 


1.567 


12.634 


484.419 


46.62 


-0.555 


X2,2o(a) 


1.000 


6.900 


375.004 


45.57 


-0.717 


X3,2o(a) 


0.957 


9.400 


447.664 


43.98 


-0.692 


X4(o) 


5.360 


2.934 


265.890 


42.89 


-0.535 


X5(a) 


5.400 


3.070 


273.490 


42.78 


-0.533 


Xe(a) 


5.843 


3.052 


273.665 


42.93 


-0.546 


Xr(a) 


6.111 


3.063 


294.382 


41.72 


-0.548 
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